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Q 1 Introduction 

U. 

A nonlinear ring cavity driven by an external signal is among the most popular models of 
nonlinear dynamics [1,2]. A special attention to this system is paid in non-linear optics where 
it was considered for the first time by Ikeda and co-authors [1-5]. Many interesting non- linear 
phenomena such as bistability, transitions to chaos, and pattern formation have been studied 
both theoretically and experimentally (see for details [6-8]). Those researches are of great 
■^f ■ interest because of prospects to apply bistable cavities as logical elements in optical computers 



and in cryptography [6,7,9,10]. We note that similar systems can be also produced not only 



Cf) \ in optics, but also in microwave diapason. 

Schematic drawing of studied device is presented in Fig. 1. A nonlinear dielectric with Kerr- 
type nonlinearity is placed inside a ring optical cavity. The cavity is excited by an input signal 
Eq (t) from an external laser source. However, we do not restrict ourselves to the nonlinear 
optical system. Our analysis is valid for any type of waves propagating in a nonlinear medium 
with cubic phase nonlinearity when output signal is partly recirculated to the input through a 
,£h ' feedback loop. 

^ ■ For the case of harmonic input signal the behavior of the system has been studied thoroughly 

[1-8,11]. Under several assumptions (plane wave approximation, neglecting dispersion and 
losses in the nonlinear media, etc.) a simple return map can be derived to describe the dynamics 
of complex amplitude of the signal 

A n+ i = A + pA n exp (i \A n \ 2 + itp) . (1) 

Here A n is the signal amplitude on n-th spreading through the cavity, A is a constant 
amplitude of the input signal, p is an amount of feedback (0 < p < 1), and <p is a total phase 
shift of the signal for a single spreading through the feedback loop. The map (JIJ is known 
as Ikeda map. For sufficiently high input power the steady-state single frequency oscillations 
become unstable and periodic pulsations arise with the period equal to doubled delay time [1-8]. 
This phenomenon is called Ikeda instability. The subsequent increase of the input signal leads 
to transition to chaos by Feigenbaum scenario and then to formation of Ikeda attractor [3,4] 
that corresponds to the fully-developed, strongly irregular chaotic motion. 

In this paper, we investigate nonlinear dynamics of this system when the input signal is 
multi-frequency, containing several discrete spectral components u>i,u>2, .-,u) n : 
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Eq = VJ A j exp (iujjt) + c.c, 



(2) 



where Aqj are complex constant amplitudes, c.c. denotes complex conjugate. For the sake 
of simplicity we suppose the frequencies Uj to be non-resonant, i.e. their ratios are irrational 
numbers, and restrict ourselves to double- frequency driving, n = 2. Note that similar problem 
(for n = 2 or 3) have been already considered in [12]. However, in cited work the main concern 
was with the steady-state characteristics of the device in the case of symmetrical pumping 
(i.e. equal input intensities for both frequencies). In contrast to [12] we will focus on case 
when the pumping is strongly asymmetric, treating the second signal as a small control action. 
Our aim is to demonstrate that varying the intensity of the second signal opens a possibility 
to effective control of the dynamics of the primary signal. We also will concentrate on the 
transitions to chaos and on peculiarities of the chaotic oscillations. Some preliminary results 
presented in [13, 14] show that the chaotic dynamics can be much more diverse than for the 
single-frequency driving. 
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Figure 1: Schematic drawing of a nonlinear ring cavity under external driving 



The paper is organized as follows. In Sec. 121 the derivation of the coupled Ikeda maps is 
described in details. Although the equations are obtained for a general case of an arbitrary 
number of input frequencies, further the double-frequency driving is considered. The stability 
analysis of obtained equations is performed in Sec. EJ Results of numerical simulations of 
non-stationary oscillations are presented in Sec. 0] Various types of transitions to chaos are 
discovered while the input signal power is increased. Sec. contains some concluding remarks. 



2 Derivation of coupled Ikeda maps 

To describe the propagation of a multi-frequency wave packet in a nonlinear medium we use the 
Hamiltonian formalism that is one of the most powerful tools in the theory on non-linear wave 
interaction (for review see [15]). We assume that only one type of waves should be considered, 
with dispersion relation u = u (k) that does not admit three-wave resonant interactions. Thus, 
first-order nontrivial nonlinear processes are caused by four-wave interactions, and Hamiltonian 
of the system can be written as follows [15]: 

H = Hq + H int + (3) 
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H = J uj (k) a,ka* k dk, (4) 

H^t = r J 7 fcifefc 3 fc4 a fc 1 a fc2 a fc3 a fc4^ {fa + h - h - fc 4 ) dk 1 dk 2 dk ?t dk 4 . (5) 

Here Hq is a linear part of Hamiltonian, Hi nt describes nonlinear interactions. In ©-© dk 
is a Fourier amplitude of the wave with wavenumber k, 5 is the Dirac delta-function, Tk 1 k 2 k 3 k4 
is a four- wave interaction matrix element (see [15] for more details). 

For the input signal (J2J) we suppose the field to be presented as a multi-frequency wave 
packet consisting of several discrete spectral components 

n 

E = A/ (x, t) exp {iujjt — ifax) + c.c, (6) 
i=i 

with frequencies ujj and wavenumbers kj being connected with each other by dispersion relation 
uj = uj (kj) and Aj being the slow varying (in comparison with complex exponents) envelope 
amplitudes. If all the frequencies are non-resonant, i.e. none of the conditions 

E^ = o 

3 

holds true for integer Nj, one can show that propagation of the wavepacket © is governed by 
the system of coupled non-linear Schrodinger equations [16] 

JdA > 'V^^+Um^wAa,^ (T) 



dt 9 dx 2 dx 2 



Here V j = % 

9 dk 



are the group velocities, uj'- = are the parameters of group 

k — \ k — k 4 



velocity dispersion, and coefficients (3 responsible for non-linear interaction between waves with 
different components of the wave packet are defined as follows: 

j ^ 'k jkjkj kj *j ^ % ^k,j k<i kj k^ I 'kj k^k^ kj I 'k^ k j k j k^ I 'k^ kj k^k * ^ ) 

Further we shall neglect the effects of group velocities dispersion (uj" = 0, V£ = V^). Hence, 
the equations (J7|) transforms into 



: ) • ( hi -h 2 ■ E b A > ') 



\ dt dx 
One can easily find the solution of © in the form 



Aj{x,t) = Aj 0,t 



x 

y 9 . 



exp 



IX 
Vn 



i \' t ~v„ 



Ai 0, t 



a; 

Vn 



• (10) 



Since the total field at the input of the non-linear medium (at x = 0) is composed by the 
external driving signal and the signal traveled through the resonator (see Fig. 1) the boundary 
condition is the following 
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E(0,t) =E + RE (L,t-t fb ). 



11) 



Here tfb is the time of signal propagation in the feedback loop, L is the length of the 
non-linear medium, and R is (generally complex) feedback parameter which is assumed to be 
the same for all frequencies. For the case of an optical ring cavity shown in Fig. 1 R is the 
reflectivity of the output mirror (other three mirrors are assumed to have R = 1). 

Substituting (J2J), ©, and (fTUj) into (fTTj) we obtain a system of delayed equations (cf. [5-8]) 



A j ( t + T ) = A oj + P A j (*) ex P 



Vn 



^(*)I 2 + £^IA(*)I 2 + ^ 



(12) 



where r = tfb + l/V g is the time delay, p — \R\, <pj — Arg (R) — lujt are total phase shifts for a 
single spreading through the feedback loop, that in general are different for different spectral 
components. Note that the phase space of the system of time delayed equations (or functional 
maps) (fT2*j) is infinite-dimensional. Our next step is to turn to a more simple system of coupled 
return maps. Examining the output signal dynamics in discrete moments of time t n = nr yields 
a system of coupled Ikeda maps 



A n + t 



A 0j + pA n , exp 



V n 



A n 



i# V 9 



(13) 



where A™ 



Aj (t r 



We restrict ourselves now to the investigation of two-frequency external signal. In that case 
we obtain from (|13|) 



A n+l 

A? +1 



Am +pA"exp 
A 02 + pA% exp 



/ln|2 PuL | I An\2 gigX , 
A l\ Vg + \ A 2\ Vg + ^ 
( I /ln|2 I3 22 L | I An\2 P21L , 



Let us suppose all four-wave interaction matrix elements to be equal: T) 



kik2ksk4 



(14) 



T. Thus, 



equations (jHJ) give /J^- = f3 = — 2T, = — 4T = 2/3. Finally, to minimize the number of 
significant parameters we make in (|T4^) the following change of variables: Aj — > Aj I \j (3L/V g , 
A ' 
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As a result we obtain 



n+l 



A n 2 +1 



A m + pA™ exp 
A 02 + pA% exp 



|^| 2 + 2|^| 2 + <p 1 
|^| 2 + 2|^| 2 + <p 2 ' 



(15) 



Since the variables A™ are complex, the phase space of this system is four-dimensional. Note 
that in the absence of the second signal, A 02 = 0, (fTTjj) can be reduced to the single Ikeda map 
©. 

We recall that coupled return maps (|TK|) derived in El performs simple model (we used 
plane wave approximation, neglected dispersion and losses in the nonlinear media, etc.). Exact 
solution of coupled Schodinger equations (JJJ) is also possible and will be presented hereafter. We 
will not discuss here the difference between two approaches either. Some remarks concerning 
differential equations with time-delayed feedback can be found in, e.g. [17-20] 
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3 Stationary regimes and stability analysis 



As stationary regimes we name henceforth the oscillation when signal amplitudes are time- 
independent. That corresponds to stable points of coupled Ikeda maps A™ = A® = const. For 



stationary intensities L 



A° 



2 



we have 



j £01 (1Q) 

1 ~ l + p 2 -2pcos$i' 1 ' 

J ° 2 (17) 



1 + p 2 - 2pcos$ 2 ' 

here I j = \A j\ 2 - intensities of outcoming signals; $j = Ii + 2Ij+ipi, i,j = 1, 2. Equations (flfij). 
(|T7jl are transcendental and can be solved only numerically. Fig. 2 represents the dependence 
h (loi) — gain characteristics. 

Dashed lines correspond to Iqi/(1 ± p) 2 , that is the limits of Ji^: 

T-^TI < h < (18) 
(1 + p) 2 3 (1-p) 2 

Fig. 2a corresponds to the case when the second signal is absent, /02 — 0, and is analogous 
to well-known gain characteristics of nonlinear ring cavity presented in, e.g., [7]. That diagram 
evidently illustrates the phenomena of multi-stability and hysteresis, peculiar to discussed sys- 
tem. The branches with negative incidence comply to unstable states. Arrows denote shifts 
from one stable branch of gain characteristics to another. 

The increase of the second signal intensity leads to the complification of I\ (ioi); number of 
branches also increases (Fig.2c). We note, that there exist stable part of gain characteristics , 
that can not be arrived by smooth varying of Joi (they are marked with dots). To obtain these 
regimes one should either provide very precise starting conditions or move in parameter space 
in very particular direction. Thus nonstationary dynamics of the system is expected to be much 
more diverse comparatively to single Ikeda map. We emphasize that in all cases performed in 
Fig. 2 the intensity of the second component remains rather small sofar it can be treated as 
weak control action (7 02 < 0.36, compare with J i)- 

Now let us undertake the stability analysis of stationary solutions. Linearizing the maps in 
the vicinity of a stable point after simple but bulky transformations we derived characteristic 
equation: 



[p 2 - 2pp (cos $i - h sin $i) + p 2 ] [p 2 - 2pp (cos $ 2 - h sin $ 2 ) + p 2 ] = 

= 16p 2 p 2 /i/2sin$isin$ 2 . 

here p — multiplicator. We note that terms in brackets in the left part of (JT§j) coincide with 
characteristic equations for single Ikeda map ((TJ. Consequently we interpret the right-hand 
term clS db coupling parameter. 

The loss of stability can happen in a various ways (see for details, e.g., [21]): as a result 
of doubling period bifurcation (p = —1), tangent bifurcation (p = 1), as well as Neimark 
bifurcation (the last corresponds to two comlex conjugate multiplicators p± = exp(±z0)). 
We plot the lines of different bifurcations on a parameter plane (A x,(fi), analyzing equation 
(fTTjj) numerically. We note that negatively inclined branches are unstable owing to tangent 
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Figure 2: Gain characteristics for coupled Ikeda maps. <fi=TT, Aq2=0 (a), Aq2=0A (b), ^4o2=0.6 (c), 
A) 2 =0.6, v?i=3vr/4 (d) 

bifurcation. Thus lines \i — 1 in a parameter plane limit the regions of multi-stability and 
hysteresis. 

Results of numerical solution ()19|) are performed in Fig.3. Solid lines denote first doubling 
bifurcations, interrupted lines correspond to tangent bifurcations, and dotted lines - to Neimark 
bifurcation. Vertical lines correspond to directions along which gain characteristics, presented 
in Fig.2, were plotted. Letters A,B,...H mark rough transitions between different branches. 
As seen from Fig.3 they coincide with hysteresis region, i.e. tangent bifurcation. Varying of the 
second signal phase </2 2 does not influence dramatically the dynamics of the system, therefore 
we restrict our numerical analysis to the case when (f 2 = 0. All results are periodic in ipi with 
period equal to 2n. 

Fig. 3a corresponds to Ikeda map. One observes the borders of hysteresis region and 
doubling bifurcation line. When amplitude of the second signal is rather small, there are no 
noticeable changes in a parameter plane. The last can be observed only when Aq2 increases to 
0.4. There appear the regions of quasi-periodic motion inside the region, limited by Neimark 
bifurcation lines (fig 3c). With the increase of /02 they grow in size, and there exist two those 
regions, dwelling in different sheets of a parameter plane (fig 3d). 
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Figure 3: Lines of different bifurcations on a parameter plane: first doubling bifurcation (solid lines), 
tangent bifurcation (dash lines), Neimark bifurcation(dotted lines). Parameters: p = 0.5, (f2 = 0, 
and Aq2 = 0.0 (a), 0.4 (b), 0.6 (c), 0.8 (d). Thin lines correspond to direction along which gain 
characteristics on Fig. 2 were built 



4 Results of Numerical Modelling and Discussion 

To investigate nonlinear dynamics of the coupled Ikeda maps f!15|) we performed numerical 
simulation in a wide range of parameters. Our primary goal is to demonstrate that adding 
the second signal component could cause qualitative changes in the dynamics of the system 
in comparison to the single-frequency driving. In Fig. 4 the typical results are presented as 
bifurcation charts (phase diagrams) in the plane (v4 i,y?i) for different values of A Q2 - The 
choice of the parameters Aqi and ip\ (that have meaning of external driving magnitude and 
frequency, respectively) is quite natural for non-autonomous nonlinear oscillator. The pictures 
are periodic in (pi with the period 2ir. The phase of the second signal (f2 does not strongly 
affect the dynamics of the system and leads mainly to shift of the bifurcation lines along the 
(fi axes. Henceforth all the presented results correspond to the case ip 2 = 0. 

In Fig. 4 the regions of periodic oscillations with different periods, quasi-periodic oscillations, 
and chaotic ones are shown. In each point on the parameter plane we determine the type of 
dynamical regime after sufficiently long transient process. For that purpose we use phase 
portraits, power spectra, and spectra of Lyapunov exponents. In the regions of periodic motion 
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all the exponents are negative, for quasi-periodic motion there are one zero and three negative 
exponents, while in chaotic regimes one of the exponents becomes positive [?,2,21]. We chose 
the direction of scanning on the parameter plane to be bottomup and from left to right and 
use the hereditary of initial conditions. The last means when the transient process is over and 
the dynamical regime in a point on the parameter plane is determined we increase smoothly 
the value of the phase shift (pi taking the current values of variables as initial conditions for a 
new point. When <pi reaches the limit (pi = 2tt we increase slightly the value of A 0i and return 
to (pi — 0. Note that owing to the bistability the bifurcation charts are multi- folded, and only 
the parts corresponding to the lower branch of the cavity transfer function are shown in Fig. 4. 
To obtain the other parts of the charts one should vary the scanning direction. 

The bifurcation chart for A Q2 = (Fig. 4(a)) coincide with that for the single Ikeda map 
(JTJ), when built in appropriate coordinates (see, for example [23], Figs. 3(b),4(b)) 1 . In Fig. 4(a) 
one can see the regions of period doublings and transition to chaos via Feigenbaum scenario. 
We show only the first two doubling bifurcations that are depicted with thin solid lines. The 
domains of higher order bifurcations are not shown because of their small sizes. In the region 
of chaotic motion there exist many domains of periodic oscillations, which have a form known 
as "crossroad area" [21,23]. The large ones (with periods 3 and 4) are shown in Fig. 4(a). The 
exit from crossroad areas occurs via Feigenbaum scenario. 

For small enough amplitude of the second signal Aq 2 (Fig. 4(b)) the changes in the bifur- 
cation chart compared with Fig. 4(a) are hardly noticeable. Only slight deformation of the 
borders of the oscillations of different types is observed. 

With increasing of A 02 the bifurcation lines become more complicated. Furthermore, in 
contrast to a single Ikeda map, quasi-periodic motion can take place (Fig. 4(c,d)). The quasi- 
periodic domains contain a large number of synchronization tongues; the largest ones are shown 
in Fig. 4(c,d,e). Inside the tongues the motion is periodic with rational winding numbers (the 
values of the winding numbers are indicated in Fig. 4). We emphasize that the winding number 
changes smoothly while moving in the quasi-periodic domains. For example, in the domain 
adjoining the region of periodic oscillations with period 1 in Fig. 2(c), the winding number varies 
from 1 to 1/2, and there exist all synchronization tongues corresponding to the rational winding 
numbers within the limits of variation. For sufficiently large A 02 (Fig. 2(d,e)) quasi-periodic 
route to chaos (Ruelle-Takens scenario) becomes possible. 

It should be pointed out that variables A 12 are the envelope amplitudes of spectral compo- 
nents with frequencies u 12 . Since a stable fixed point of the map ([T5]) describes a steady-state 
solution when the amplitudes are time-independent, Ai >2 = const, the real physical field E (x, t) 
(jHJ) in general case (for the frequencies u>i >2 supposed to be non-resonant, see Sec. |2J) exhibits 
quasi-periodic oscillations with two independent frequencies, in contrast to the single Ikeda map 

where steady-state solution corresponds to periodic oscillations with frequency of external 
driving. In non- stationary regimes when the amplitudes Ai j2 vary in time, both components of 
the output signal become modulated. We presume the carrier frequencies uj\ )2 to be sufficiently 
far off from each other {\u)\ — uj 2 \ ^> 2n/r, r is the time delay), so their modulation spectra do 
not intersect. So, dealing with an experimental device, it is possible to investigate the dynamics 
of each component separately using a band-pass filter tuned on one of the carrier frequencies. 

When the first period doubling bifurcation occurs (appearance of a 2-period cycle from a 

x In Ref. [23] the bifurcation charts are built in the plane of parameters that correspond to A and 1/ip in 
the equation Q). After appropriate redrawing one can demonstrate that the maps become 27r-periodic in if. 
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Figure 4: Bifurcation charts on the plane of parameters (Aq±, (pi) for p = 0.5, if2 = 0, and A02 = 0.0 
(a), 0.2 (b), 0.6 (c), 0.8 (d). Transitions to chaos via Feigenbaum scenario are shown by dashed 
lines; quasi-periodic transitions to chaos are shown by dotted lines; period doubling bifurcations are 
shown by thin solid lines Only first two bifurcations are depicted, higher-order ones are not shown 
for the lack of space. Ch — regions of chaotic oscillations Q — regions of quasi-periodic oscillations, 
numbers mark regions of periodic oscillations with corresponding period. For synchronization tongues 
the values of winding number are indicated. 

stable point), the amplitudes exhibits periodic self-pulsations with period 2r. Therefore, this 
bifurcation corresponds to the oscillatory Ikeda instability caused by the delayed feedback [1-5]. 
The spectrum contains side-bands (satellites) with frequencies uip ± nu>i, n = 1, 2, . . ., where 
uji = tt/t — the self-modulation frequency. During transition to chaos the spectrum is enriched 
by new sub-harmonic components with frequencies uoi^ ± nuoi/2 m . 

On the contrary, in the quasi-periodic regime a fixed point loses its stability via the Hopf bi- 
furcation, when a unit circle is crossed by a pair of complex conjugate multipliers and an attrac- 
tor is invariant curve [1,2,21]. The side-band frequencies are lo\^ ^nujq with ujq ^ n/mr. For 
the dynamics of the total field E (x, t) (JUJ) it corresponds to periodic amplitude self-pulsations 
with period different from 2r, and the attractor is a 3D torus with three independent frequen- 
cies (two career frequencies Ui t2 and that of the self-pulsations, ujq). Inside the synchronization 
tongues ujq is commensurable to the Ikeda frequency ujj, i.e. a winding number w = ujq/uj 
is rational, w = p/q, p and q — integer. The single Ikeda map (pQ) does not exhibit the Hopf 
bifurcation and quasi-periodic motion [1-5,8,21,23]. Hence, we obtain a new type of insta- 
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bility caused by nonlinear interaction of two driving signal components, i.e. by phase cross- 
modulation. Note that the effects of phase cross-modulation are very important in nonlinear 
optics, particularly for wave propagation in optical fibers [24] . 

We point out that the results of direct computer modelling of (0) are of great agreement with 
the analysis of characteristic equations undertaken in Sec. 01 Some distinctions are connected 
with the choice of scanning direction. 

We observe three different routes to chaos: Feigenbaum period doubling scenario, quasi- 
periodic (Ruelle-Takens) scenario, and hard transition. Accordingly, three types of chaotic 
attractors exist (Fig. 5). We use the following notation: Cp for the Feigenbaum attractor 
originated from period doubling bifurcation sequence, Cq for the one formed from quasi-periodic 
motion, and Cj for the well-known Ikeda attractor [5] that exists for high pumping levels and 
corresponds to fully-developed, strongly irregular chaotic motion. The three described types 
differ not only in visual appearance, but also in signature of Lyapunov exponents: Cg-type 
attractor has one positive, one zero, and two negative exponents, while Cp and Cj types 
have one positive and three negative exponents. No hyperchaotic regimes with two positive 
exponents occur. We emphasize that for the single Ikeda map only Cp and Cj types exist 
and quasi-periodic route to chaos never takes place. We represent in Fig. 5 Feigenbaum-type 
attractor (a), Cp (<pi = it/6, A m = 1.3); (b) attractor originated from the quasi-periodic route 
to chaos, Cq (<fi = 3n/2, A m = 1.5); (c) Ikeda-type attractor, Cj (ipx = it, A m = 2.5). Other 
parameters are the same as for Fig. 4(c). 




1.0 



0.5 - 



0.0 



-0.5 
0.0 





2.0 








1.0 




0.0 




-1.0 




-2.0 




(c) 



0.5 1.0 1.5 2.0 0.0 1.0 2.0 3.0 4.0 5.0 

Re (A n ) 



Figure 5: Chaotic attractors of different types 



From Fig. 3 one can see that the bifurcation charts are very complicated, and thus, a 
great variety of bifurcation sequences could be observed while moving in different ways on 
the parameter plane. For better understanding the behavior of the system we built single- 
parameter bifurcation diagrams: Re (A^) plots versus Aq\. Some examples are presented in 
Fig. 6 for different values of <pi. All diagrams are plotted with smooth increase of A i with 
hereditary of initial conditions for (pi = tt/Q (a), (fi = n (b), and (pi = 3ti/2 (c). Other 
parameters are the same as for Fig. 4(c). 

The first diagram (Fig. 6(a)) shows period doubling bifurcations and transition to chaos 
via Fiegenbaum scenario at A m pa 1.2. In the region of chaotic motion several windows of 
periodicity are clearly seen. At A m pa 1.5 the Cp-type attractor begins to transform in the 
Cj-type one. When A$\ pa 1.63 one observes a hard transition between different types of 
chaotic oscillations known as a crisis of chaotic attractors. In that case chaotic attractor 
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Figure 6: Single-parameter bifurcation diagrams -Re versus Aqi 



abruptly expands or contrarily shrinks, that is explained by its collision with non-stable periodic 
orbit [21,22]. Further increase of Aqi leads to the reverse doubling bifurcation sequence. At 
Aqi w 1.8 a hard transition to quasi-periodic motion takes place. Inside the region of quasi- 
periodicity one can see several windows of synchronization. With increase of Aqi the quasi- 
periodic attractor shrinks, and at Aqi « 2.0 the motion again becomes periodic. Finally, the 
hard transition to Ikeda attractor Cj takes place at A i « 2.25. 

On the second diagram (Fig. 6(b)) one can see a region of quasi-periodic oscillations that 
appears and disappears softly (i.e. via direct and reverse Hopf bifurcations) at Aqi 0.8 
and Aqi ~ 1.1, respectively. Hard transition to chaotic attractor Cj occurs at Aqi ~ 1.25. 
Subsequent increase of Aqi causes the expansion of the attractor Cj. 

The bifurcation diagram presented on Fig. 4(c) illustrates soft appearance of quasi-periodic 
oscillations (Aqi ~ 0.9). A number of synchronization windows with different periods can be 
observed, inside some of them period doublings takes place. At Aqi ~ 2.0 a crisis of chaotic 
attractors is clearly visible. The bifurcation sequence ends with a hard transition to Ikeda 
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attractor at Aqi ~ 2.3. 

5 Concluding remarks 

In this paper, we discussed an extension of the Ikeda system for the case of multi-frequency 
external signal. The system of coupled Ikeda maps was derived to describe the dynamics of 
slowly varying amplitudes of spectral components. For the simplest case of double-frequency 
driving results of numerical simulation were presented, that confirm that the dynamics is much 
more diverse than for the single-frequency driving. Various transitions to chaos were observed: 
the Feigenbaum period doubling scenario, quasi-periodic (Ruelle-Takens) scenario, hard tran- 
sition to chaos, and crisis of chaotic attractors. Applying of the second input signal gives the 
opportunity to efficient control of the dynamics of the system. Changing amplitude and phase 
of the second signal component one can either suppress the chaotic oscillations or conversely 
cause the steady-state oscillations to be unstable. These results are significant for prospects to 
use nonlinear ring cavities as logical elements in optical computers [7]. 

However, the validity of the model (j!5|) needs further consideration, since a number of 
simplifying assumptions has been made. First of all, we neglected transverse effects (diffraction, 
self-focusing, etc.) that lead to a wide class of transversal instabilities [7,8]. We also did not 
take into account group velocity dispersion effects that may cause modulation (Benjamin-Feir) 
instability and formation of solitons [7,24]. All these phenomena usually play an important 
role in nonlinear optics, and their consideration will be a subject of further investigation. 

Finally, we note that Ikeda map describes non only a nonlinear ring cavity, but also many 
other interesting physical systems, e.g. a periodically kicked nonlinear oscillator [23]. Evidently, 
the system of coupled maps (fTKjl can be derived for a quasi-periodically kicked oscillator. 
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